Long-time Protection of Nonlocal Entanglement 
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We investigate how nonlocal entanglement, as identified by violations of a Bell inequality, may 
be protected during the evolution. Our system consists of two qubits each embedded in a bosonic 
reservoir evolving independently and initially in an entangled mixed state. We show that the 
violation of the Bell inequality can be related to the single-qubit population of excited state in such 
a way that, by appropriately choosing structured environments that give rise to sufficiently high 
values of population trapping, long-time protection of nonlocal entanglement can be correspondingly 
achieved. 



I. INTRODUCTION 

Entanglement is one of the fundamental but coun- 
terintuitive aspects of quantum mechanics [3| ■ A ba- 
sic trait of entangled, or quantum correlated, systems 
is that a measurement on one subsystem of the com- 
posite system determines the state of the other, inde- 
pendently on their separation. Entangled systems can 
also be characterized by the presence of quantum corre- 
lations that cannot be simulated by any classical local 
modelf^]. However, while for pure states the presence 
of entanglement always corresponds to the presence of 
nonlocal quantum correlations (NLQCs)[l, 01, this is not 
the case for mixed states; in fact it has been demon- 
strated that there are bipartite entangled mixed states 
whose correlations admit a classical local model [5] . These 
NLQCs are unambiguously identified by violations of Bell 
inequalities 0, Q. These violations have been confirmed 
by several experiments 0, H, and are crucial for some 
applications in quantum information (lo| . such as to guar- 
antee the safety of device-independent key distribution 
protocols in quantum cryptography [Til IT3| . 

On the other hand, realistic quantum systems present 
unavoidable interactions with their environment [l3[ that, 
during the evolution, give rise both to decoherence and 
to entanglement losses, turning any initial pure state into 
a mixed one. Although these two processes appear to 
be related, it has been shown that two initially entan- 
gled qubits, each locally interacting with a correspond- 
ing environment, formed by a zero temperature memo- 
ryless (Markovian) bosonic reservoir, may be subjected 
to "sudden death" In this process, the qubits be- 
come disentangled at a finite time, in contrast to the 
single qubit exponential decoherence. This phe nomenon 
has been recently experimentally confirmed[l5l [l6| . In- 
terest therefore arises in finding physical conditions and 
methods able to protect entanglement. Recently, meth- 
ods aimed at extending the entanglement storage time 
have been proposed. They make use of environments 
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with memor y ( non-Mar kovian) , where entanglement re- 
vivals occur [TtI ITsI I. of quantum Zeno effect . or of 
structured environments like photonic crystals, where en- 
tanglement can be trapped [20l|. 

We shall consider a system composed by a couple of 
qubits, each embedded in a bosonic reservoir, evolving 
independently to examine the connection between viola- 
tion of a Bell inequality and single-qubit excited state 
population. We shall study the special case when the en- 
vironment is such that sufficiently high population trap- 
ping occurs (as in photonic-band gap materials 21]), to 
see if a long-time protection of nonlocal entanglement, as 
identified by Bell inequality violations, is achievable. 



II. FORMALISM 

We begin by considering a bipartite system made of 
two identical independent two-level atoms (qubits) each 
embedded in a zero temperature bosonic environment. 
We indicate the two qubits with A and B, each having 
ground and excited states |0), |1) separated by a tran- 
sition frequency loq. A crucial point for our study is to 
obtain the two-qubit reduced density matrix at the time 
t. This can be accomplished by exploiting a procedure 
based on the knowledge of the single-qubit dynamics [13] ■ 
The fact that the two parts "qubit -|-environment" evolve 
independently allows to write the time evolution oper- 
ator of the composite system as a tensor product of 
two terms relative to the two parts. This in turn leads 
to the fact that, given the time-dependent single-qubit 
density matrix elements as pf^,{t) = XliZ' ^fi' (0), 
pfA^) = T.mm'B]]V'\t)pl^,{Q), the time-dependent 
two-qubit density matrix elements result of the form[l7j 
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where the indexes i, j,l,m — The condition of inde- 
pendent evolution is applicable also to the case of com- 
mon environment provided that the atoms are at a dis- 
tance larger than the spatial correlation length of the 
reservoir and such that the dipole-dipole interaction can 
be neglected. 
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For each part, the dynamics of the qubit-environment 
interaction is described by the Hamiltonian [T3| 

H = hujoa+a^ + ^ ^hhJkblbk + (^gk^+bk + gla^b^^ , 

k 

(2) 

where a+ = |1)(0|, u- ~ |0)(1| are the qubit raising and 
lowering operators, bj., bk are the photon creation and 
annihilation operators and gk is the coupling constant of 
the mode k with frequency ojk- 

When the environment is at zero temperature the 
single-qubit reduced density matrix p^{t) can be writ- 
ten, in the basis {|1), |0)}, asfisf 

V PlMq*{t) pfo(0)+pfi(0)(l-|g(t)n/ 

(3) 

where p{t) = pfi(0)|g(t)P represents the single-qubit 
population of excited state at time t and is the 

population parameter. The previous equation shows that 
the single-qubit dynamics depends only on q{t). The two- 
qubit density matrix p(t) is then obtained by means of 
Eq. JT]) , its elements depending only on their initial values 
and on the function q{t). In the following these density 
matrix elements will be meant in the standard two-qubit 
basis B = {|1) = 111), |2) = |10), |3) = |01), |4) = |00)}. 

III. CHSH-BELL INEQUALITY 

In order to individuate when the entanglement is non- 
local we have to choose a Bell inequality and look for 
its violations. Here we shall exploit the well-known 
Clauser-Horne-Shimony-Holt (CHSH) form of the Bell 
inequality^]. The CHSH-Bell inequality is the most ex- 
tensively used among the Bell inequalities. The rea- 
son of this relies on the fact that it is experimentally 
testablej6i, l^] and that for the case of two spin-like 
measurements (two measurements with two outcomes 
per each) any other Bell inequality is equivalent to the 
CHSH-BeU one^]. 

Let the operator Os = Os{0s,(j>s) be a spin-like ob- 
servable with eigenvalues ±1 associated to the qubit 
S — A,B, defined as Os — O5 • 0-3, where Og = 
(sin 6*5 cos (/)s, sin 6*5 sinew's, cos 6*5 ) is the unit vector in- 
dicating a direction in the spin-like space and 0-3 = 
(erf, erf, af) the Pauli matrices vector. The expression 
of Os in terms of the basis states is 

O3 = cos05(|l)(l|-|O)(O|) 

+ sin05(e''^-11)(O|+e"'^-lO)(l|). (4) 

The CHSH-Bell inequality associated to the two-qubit 
state p for the operator O is@ 

Bip) = \{OaOb) - {OaO's)\ + {O'aOb) + {O'M < 2, 

(5) 



where B{p) is the Bell function, {OaOb} = TtIpOaOb} 
is the correlation function and O'g = 03{9'g,<j>s). If, 
given the state p, it is possible to find a set of angles 
and {eA,0'A,eB,0'B} such that the CHSH-Bell 
inequality is violated, that is B{p) > 2, the correlations 
are nonlocal and cannot be reproduced by any classical 
local model. 

The Bell function at time t is then obtained as a func- 
tion of q{t) from the evolved two-qubit state p{t). Our 
aim is to individuate CHSH-Bell inequality violation re- 
gions characterized by B{p) > 2. In particular, we con- 
sider the case when the two-qubit density matrix at time 
t has an "X" structure, that is having non-zero elements 
only along the main diagonal and anti-diagonal. This 
has been shown to happen if the initial two-qubit density 
matrix has itself an "X" structure: in fact this form is 
preserved during the time evolution due to the Hamilto- 
nian H — Ha + Hb, where Hs {S = A,B) is the sin- 
gle qubit-environment Hamiltonian of Eq. ([2]) [17]. The 
choice of such initial conditions relies on the fact that, on 
one hand, it simply leads to analytical solutions and, on 
the other hand, it includes several initial states of phys- 
ical inter est, as Bell-like (pure) or Werner-like (mixed) 
states [l0l.[23j. 

The maximum of the Bell function B{p) of Eq. ^ 
is found using the standard procedure 3] at the angles 
{01 ^2} = {{k + k')TT - (5i4 - (^23, {k - k')TT - Si4 + S23}, 
where k,k' are integer numbers and (5i4, S23 are re- 
spectively the phases of pi4(0), ^23(0) and in addition 
{0i,e[,02{t),e'^{t)} = {0,f,arctan ^,^-02(01, and 
in this case results to bef23] 

i?max(/5(i)) = 2^V^t) + Q^{t), (6) 

where 

V{t) = 1 - 2\qit)\^[l + pn{0) - P44(0) - 2pn(0)k(i)P], 
Q(t) = 2|g(<)p(|pi4(0)| + |p23(0)|). (7) 

These equations immediately highlight the relation be- 
tween Bell function and population parameter |(?(i)p. We 
notice that this expression for -Binax(/3(i)) coincides with 
that given by the formal Horodecki criterion [2^. 

IV. STUDY OF CHSH-BELL INEQUALITY 
VIOLATIONS 

In this section we study for which values of |(7(t)P 
CHSH-Bell inequality violations occur. 

A. Initial states 

We consider as two-qubit initial states extended 
Werner-like (EWL) stat esjl8| that are mixed states re- 
ducing to Werner states Q or to Bell-like states under 
certain conditions. In this way we shall be able to analyze 
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the effects of both mixedness and degree of entanglement 
of the initial states on the dynamics of the maximum 
of Bell function -Bmax(/5(i))- Moreover, Werner states 
are encountered during experimental generations of Bell 
states S3 ■ EWL states are defined as 



/5*(0) 
/5*(0) 
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r|1')(*| 
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(8) 



where r indicates the purity of the initial states, I4 is the 
4x4 identity matrix and 

1$) =Q|01) + ^e'*|10), I*) =a|00) + ^e'*|ll), (9) 

are the Bell- like states with a, (3 real and + = 1. 
The EWL states of Eq. ^ reduce to the Werner-like 
states for a = ±/3 = I/a/Z, that is when their pure part 
becomes a Bell state. For r = the EWL states become 
totally mixed states, while for r — I they reduce respec- 
tively to the Bell-like (pure) states |$),|*) of Eq. 
The EWL states /5*(0), /5*(0) are more general than the 
Werner-like states since they can have a non maximally 
entangled pure part. 

The initial two-qubit density matrix elements for the 
state /5* are 



pfi(O) 



1 



P22(0) 



1 
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1-^ 2 
— ^ 1- a 



l-r 



0, p* (0) = a/3e*V, 
while for the state /5* 
l-r 



(10) 



pfi(O) = 
P3*3(0) = 



4 

l-r 



ptM 



pU^) 

1 - 



4 ' ' 4 

(0) = a/3e*V, P2*3(0) = 0. 



-I- a^r. 



(11) 



These density matrix elements of the EWL states are 
such that the resulting density matrix has an "X" struc- 
ture. 



B. Nonlocal entanglement protection 

In order to evidence the regions of violations of the 
CHSH-Bell inequality, we plot the maximums of Bell 
function -Bmax(p*(0) ^'i'^ S,iiax(p*(i)), relating to the 
two EWL states of Eq. (I8|), as functions of the popula- 
tion parameter |9(i)P for different values of the purity r 
and for fixed a. Fig. [1] relates to the EWL initial state 
/5*(0) while Fig.[2]to the EWL initial state p*(0). Both 
figures show that violations of CHSH-Bell inequality can 
be achieved {B > 2), for r assigned, when the population 
parameter |q(i)P is larger than a corresponding threshold 



Bmax(P4,(t)) 
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FIG. 1: Maximum of Bell function -Bmax{p4. (*)) as a function of 
the population parameter |gp starting from the initial EWL state 
/3*(0) with a = 13 = l/\/2 for different values of purity r: r = 1 
(solid curve), r = 0.9 (dotted curve), r = 0.8 (long-short-dashed 
curve), r = 0.7 (long-short-short-dashed curve), r = 0.6 (long- 
dashed curve). 
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FIG. 2: Maximum of Bell function -Bmax(/3a<(t)) as a function of 
the population parameter |gp starting from the initial EWL state 
p*(0) with o = /3 = I/V2 for different values of purity r: r = 1 
(solid curve), r = 0.9 (dotted curve), r = 0.8 (long-short-dashed 
curve), r = 0.7 (long-short-short-dashed curve), r = 0.6 (long- 
dashed curve). 



value |(7r-P- In particular, for initial maximally entangled 
pure states (r = 1), the threshold values are — 0.8 
for the Bell state |$) = (|01) ± |10)/V2 and jqip « 0.77 
for the Bell state 1^-) = (|00) ± |11)/V^. However, from 
figures 1 and 2 one sees that for small enough values of 
purity r (r < 0.71 for both initial states), CHSH-BeU 
inequality violations never occur. Therefore, for initial 
conditions such that the purity of the states is sufhciently 
high, if physical conditions of the environment are such 
that high values of single-qubit excited state population 
can be achieved, nonlocal entanglement between the two 
qubits is obtained. As a consequence, when these high 
values of population can be trapped, a long-time protec- 
tion of nonlocal entanglement follows. These required 
values of population trapping may be reached, for exam- 
ple, in the case when the environment is constituted by 
a photonic crystal exhibiting photonic-band gap(20l. [2]|. 
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V. CONCLUSIONS 

In this paper, we have considered a system composed 
by a couple of initially entangled qubits, each embedded 
in a zero temperature bosonic environment, that evolve 
independently. In this case, the two-qubit density matrix 
at time t is simply expressed in terms of the coefficients 
that determine the evolution of each single independent 
qubit and on the initial conditions of the qubit pair. 

In order to individuate when the two-qubit entangle- 
ment presents nonlocal correlations we have exploited 
the CHSH-Bell inequality and found its relation with 
the single-qubit excited state population represented by 
a population parameter. We have graphically shown this 
relation taking as two-qubit initial conditions Werner-like 
mixed states. We have found that, when the initial purity 
of the states is high enough, for some range of the pop- 
ulation parameter values, a violation of the CHSH-Bell 



inequality is achieved. We have suggested that, in the 
case of structured environments engineered so to present 
single-qubit population trapping, as for example photonic 
crystals presenting photonic-band gaps, long-time entan- 
glement protection can be reached. 

Nonlocal entanglement means that a Bell inequality is 
violated and that the two-qubit correlations cannot be 
reproduced by any classical local model. This condition 
appears to be relevant both for an applicative point of 
view, e.g. for quantum information purposes pJll ll^ . and 
for a fundamental point of view. Our results demonstrate 
that purely quantum nonlocal correlations can be main- 
tained during the evolution and highlight the importance 
of using appropriately engineered environments. 
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